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Abstract 

We derive the non-linear integral equations determining the free energy of the 
three-state pure bosonic Uimin-Sutherland model. In order to find a complete set 
of auxiliary functions, the anti-symmetric fusion procedure is utilized. We solve the 
non-linear integral equations numerically and see that the low-temperature behavior 
coincides with that predicted by conformal field theory. The magnetization and mag- 
netic susceptibility are also calculated by means of the non-linear integral equation. 
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1 Introduction 



The Bethe ansatz (BA) method [jlj is a most fundamental approach to the study of one 
and two-dimensional exactly solvable lattice models. In the BA strategy, the quantum 
transfer matrix (QTM) method has been found quite a powerful tool to investigate the 
thermodynamics of several spin models and electron systems 0-0. In most cases, after 
applying the BA it proved useful to introduce auxiliary functions which are ratios of the 
components of the eigenvalue of the QTM and their non-linear integral equation (NLIE). 
The main merit to utilize the NLIE is that we do not have to use the string hypothesis fi~0|| , 
which sometimes raises fundamental questions about its validity as well as pragmatical 
questions about the accuracy of the unavoidable truncation procedures of the typically 
infinitely many integral equations in the traditional thermodynamical Bethe ansatz (TBA) 

©• 

Once discovering a good (finite) set of auxiliary functions, we can get as good an 
accuracy as we need, because no truncation is necessary. Several remarkable studies of 
the thermodynamics [6|, Q and the mathematical aspects |8| have been achieved by this 
scheme. 

In this paper, we consider the one-dimensional Uimin- Sutherland (US) model [O]. For 
m fermionic and n bosonic components the system is referred to as the (m, n) model. 
For this general case the quantum transfer matrix is constructed and its eigenvalue equa- 
tions are derived. The (3, 0)-US model is studied in detail for which the spins attached 
to the sites take three states, which can be regarded as the three vectors of the vector 
representation of the Lie algebra si (3). We have succeeded in finding the NLIE for the 
(3, 0)-US model with regular analyticity properties. To obtain these well-posed NLIE for 
the (3, 0)-US model, it is crucial to consider the anti-symmetrical fusion (ASF) model in 
addition to the defining one. The ASF of the vector representation of si (3) is another 
three-dimensional one, and thus called conjugate. The spinon like picture in the (3, 0)-US 
model is made clear with the help of the NLIE. As an application some physical quantities 
like specific heats are calculated with good accuracy by numerical treatments of the NLIE. 
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The remainder of this paper is organized as follows. We give an explanation of the 
US model and the QTM method in Sec. 2. This section is also intended to fix notations. 
In Sec. 3, we apply the BA to the QTM of the general US model and we also introduce 
the NLIE for the (3,0) model. In several limiting cases, for example the SU(2)-\imit, 
exact analytic calculations are performed. In Sec. 4, we solve the NLIE numerically and 
obtain the entropy, specific heat, magnetization, and magnetic susceptibility. Some low- 
temperature properties exposed by these data are discussed. In Sec. 5 we give a summary 
of our work. 

2 The Uimin-Sutherland model and its QTM 

Let us begin with the review and definition of the general US models. The general one- 
dimensional g-state US model is defined as follows. Consider a one-dimensional lattice 
with L sites and periodic boundary conditions imposed. A g-state spin variable on is 
assigned to each site i. We can generally consider the situation where each spin a has its 
own grading, i.e. statistics number e a = (— 1)^ Q = ±1. A spin a with e a = +1 (e a = — 1) 
is called bosonic (fermionic). The Hamiltonian of the US model can be introduced as 

L 

Wo = ^7T i)i+ l (1) 
i=l 

with the permutation operator vr^j+i 

7Tj.j + i|ai • • • aj«i+i • • • «l) = a i+1 ai ■ ■ • 

where £, ai a i+1 is 1 if both «j and oti + \ are fermionic, and otherwise. 

Model (||) is shown to be exactly solvable on the basis of the Yang-Baxter equation. 
Many well-known exactly solvable models are of type ([!]), e.g. the spin-1/2 Heisenberg 
chain corresponds to q = 2 and ei = £2 = +L the free fermion model to q = 2 and 
ei = —62 = +1, the supersymmetric t-Jmodel to q = 3 and ei = —62 = £3 = +L If m of q 
e's are +1 and n(= q — m) are —1, for example, e± = ■ ■ ■ = e m = +1, e m+ i = • • • e q = — 1, 
we call the model (m, n)-US model. 
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Now we want to consider the (3, 0)-US model (q = 3, t\ = e<i = £3 = +1) whose 
Hamiltonian is equivalent to that of the SU(2) spin-1 chain defined by 

L 
i=l 

For the study of the thermodynamics we introduce the QTM as follows. First, the clas- 
sical counterpart to (|l]) the Perk-Schultz (PS) model [jl2] is defined on a two-dimensional 
square lattice of L x N sites, where we impose periodic boundary conditions throughout 
this paper. We assume that variables taking on values 1, 2, • • • , q are assigned to the bonds 
of the lattice. The Boltzmann weight associated with a local vertex configuration a, f3, /j, 
and v is denoted by R^(v), where v is the spectral parameter (Fig.l). 

Using the Yang-Baxter equation, a lattice model defined by 

Kp( v ) = S auS„ + v • (-I)^m • 6 a/3 6^ (2) 
is proved to be exactly solvable. Defining the row-to-row transfer matrix 

the partition function is given by 7-l,n = TrT N (v) where the trace is taken in the q - 
dimensional space. In this paper we are not primarily interested in the PS-model itself, 



but rather in its Hamiltonian limit obtained from T. Making use of Baxter's formula [13] 
at v = 



d 



H = —\nT{v) 

av 



22^i,i+i, (3) 



v=0 i=1 

we get the Hamiltonian of the g-state US model. 

The main idea of the quantum transfer matrix (QTM) method at finite temperature 
is as simple as follows (for details the reader is referred to the papers H). First, let us 
define a new set of vertex weights R(v ) by rotating R(v) by 90 degrees as 

and consider the transfer-matrix T which is the product of R(v). Upon introducing a 
large integer N (Trotter number), we get the following relation with arbitrary reciprocal 



temperature (3 



T(-(3/N)T(-(3/N) = e - 2 4n+o( W Nf) 



Finally, the partition function of the one-dimensional Hamiltonian TL at finite temperature 
T = 1/(3 can be calculated by means of the following "Trotter-Suzuki" formula 

z = Tve -m = lim (TrT(u)T(u)) N/2 , u = -(3/N. (4) 

N^oo 

In other words, the finite-temperature partition function for the one-dimensional Hamil- 
tonian is calculated as that of a staggered two-dimensional vertex model. For technical 
reasons, it is convenient to define another vertex weight R(v ) as a rotation of R(v) by -90 
degrees 

With these preparations, the quantum transfer matrix Tqtm corresponding to the contri- 
bution of columns to the partition function is given by 

N/2 



( T QTM^ {v) = E n R ^-_ T& _ Siv + u)R^(iv - „). (5) 



[v) = 

(i 3=1 

Here we have introduced a spectral parameter v such that 'Tqtm(^) is a commuting family 
of matrices generated by v. This will allow us to diagonalize Tqtm- The final results, 
of course, are physically interesting only for v = as the partition function of the one- 
dimensional US model at temperature l/{3 is given by 

2= lim TrT Q L TM (0). 

The free energy per unit length is 

/ = - lim In Z = -\ In A max (0), 

where A max (-u) is the largest eigenvalue of the QTM. If we succeed in obtaining the next- 
largest eigenvalue Ai(t>), the correlation length £ at the finite temperature T = 1/(3 is 
given by 



£ 1 = — lim In 



Ai 



Amax 
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3 Nonlinear Integral Equation for the Free Energy 

Let us apply the BA for this QTM. While keeping solvability, we can add external field 
terms TL ex t to TLq like 



L q 



H = Hq + Ti cxt = Tip - /I. 



The BA for the eigenvalues of the QTM of the general (/-state PS model is conjectured to 
take the form of 

A( V ) = 

M J~ 1 v-vl' 1 -ien ™i v-v{_+ie j 



— i — J- fcj 2 kj — 1 kj 



M q -i „, „,g-l 



v — vi — ie n 



{j = 2,---,q-l), 

\ q (v) = [] 1 """'l-i — ( v + iu + ( v ~ iu)^e^ q - (6) 

We concentrate on the (3, 0)-US model. Defining qi(v) = ]J(v — and </>±(«) = (t> ± 
iu) N / 2 , the last formulas are reduced to 

qiw 

qi(u -i) q 2 (u + i) , , , , . «„ 2 

A 2 («) = -—(/>+ (w)^-(v)e W2 , 

qi(f) q2(u) 

A 3 («) = q2{V ~ l) Mv + i)<P-(v)e^, 

q2W 

A( W ) = A 1 ( V )+A 2 ( V ) + A 3 W. (7) 

We have checked numerically for small iV's that the above BA gives the correct eigenvalues. 

Next, we want to define auxiliary functions which determine the (largest) eigenvalue 
completely and satisfy a closed system of integral equations. In the case of the t- Jmodel, 
a complete set of auxiliary functions is found || to be 

t Ai 7 A3 A1A3 , . 

-> b= \ — r-^' c= v-7^ — ; — \ — 1 — TT' ( 8 ) 



A 2 + A3 Ai + A 2 ' A 2 (Ai + A 2 + A3) 



Unfortunately, for the (3, 0)-US model the above set does not admit a closed set of func- 
tional equations to determine the eigenvalue. In order to complete this set we found in 
the fusion procedure [|14|] a useful working basis |7], [9||. The fusion model is defined as 
follows. The PS model given by (Q) is identified with the vector representation of the 
sl(m\n) model. On the basis of tensor products and proper projections, we get the other 
representations with higher levels. The solvability of the fusion model is again guaranteed 
by the Yang-Baxter relation. Here we consider only the ASF model. For the case of sl(n), 
we call the (n — l)-th ASF model the conjugate model because the ASF of the fundamental 
model and the (n — l)-th ASF is zero-dimensional. This conjugacy is somewhat like that 
between quark and anti-quark representations of SU (3). For the si (2) case, the ASF of two 
vector representations (i.e spin- 1/2) gives the zero-dimensional representation (i.e spin-0). 
Therefore, we call the sl(2) PS model self-conjugate. Under the conjugate transformation, 
we see 



DO 



For the si (3) case, 



m 



2 








1 








1 






2 


<— > 




3 


<— > 


3 


1 






3 
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The BA equation for the fusion model is obtained by the simple replacement 



Xi(v) -> A 2 ,3(v), X 2 (v) -> X lj3 (v), X 3 (v) -> Xx t2 (v), 
A(v)^A{v) = X 2 ,3(v)+X 1:3 (y) + X 1>2 (y), 



0) 
(10) 



where Xi >m (v) is defined by 



Xi, m (v) = X[(v)X m (v + i). 



We are ready to define the auxiliary functions. In the case of the sl{2) PS model, 
the ASF model is nothing but the fundamental model itself. As is known, a complete 
set of auxiliary functions consists of p = X\/X 2 and p = p^ 1 = A2/A1. In fact, p and p 
are conjugates. For the sl{3) case, the three functions 6, b and c are not complete in the 
(3, 0)-US model, however the 6 functions b, 6, c and their conjugates constitute a complete 
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set. For this reason we have introduced the conjugate transformation. Let us proceed this 
program. In dependence on the real variable x, we define the functions 



si(x) = 

S3(x) = 



Ai 



A2 + A3 
A3 



A12A 



12^23 



S5(x) 



Ai + A 2 
A1A3 



v=x+i/2 
v=x—i/2 



Al3(Al2 + A23 + A13) 

A12 



A 2 (Ai + A 2 + A 3 ) 
A(x) = Ai + A 2 + A 3 | w= 



■-x ' 



S2(x) 

s 4 (x) 

Al2 + A13 

A(s) = A12 + A 23 + Ms\ v=x - i/2 • 



v=x—i/2 



A13 + A23 
A23 



(11) 



After a lengthy calculation using the Fourier transform the following non-linear integral 
equation are proved explicitly 



/ In si (a;) \ 




/ f3e 1 (x) 




1 


K 


-K x -K x 


-K 3 


-K 3 


-K 4 \ 




/ In 5x (X) 


\ 


lns 2 (x) 




(3e 2 (x) 






-K 2 


K -K x 


-K 3 


-K 6 


-K 3 




ln5 2 (x) 




lns 3 (x) 




(3e 3 (x) 






-K 2 


-K 2 K 


-K 5 


-K 3 


-K 3 




ln5 3 (x) 




In S4(x) 




/3e 4 (x) 




V 


-K 3 


-K 3 -K 4 


K 


-K x 


-K x 


* 


11154(0;) 




lns 5 (x) 




/3e 5 (x) 






-K 3 


-K 6 -K 3 


-K 2 


Ko 


-Kx 




ln5 5 (x) 




V lns 6 (x) / 




V (3e e (x) 


J 


V 


-K 5 


-K 3 -K 3 


-K 2 


-K 2 


Ko ) 




\ ln5 6 (x) 


/ 



Several remarks are in order. First, Si = 1 + Si (i = 1, ■ 
terms/bare energies of the spinous are defined as 



(12) 

,6). Second, the driving 



ei(x) 


= Vx{x) 


+ (-2/xi + H2 + M3 )/3, 


(13) 


e 2 (x) 


= Vl(x) 


+ {m - 2^2 + M3)/3, 


(14) 


e 3{x) 


= Vl(x) 


+ (Mi +M2 - 2// 3 )/3, 


(15) 


e 4 (x) 


= V 2 (x) 


+ ("Mi - M2 + 2/x 3 )/3, 


(16) 


e 5{x) 


= V 2 (x) 


+ (-/xi + 2/x 2 - /U 3 )/3, 


(17) 




= V 2 (x) 


+ (2/xi - /i 2 - M 3 )/3 


(18) 



with 



VAx) 



2tt 



^/3 2cosh(2vrx/3) - 1' 
Lastly, the kernels Ki(x)(l = 0, • • • , 6) are 



V 2 (x) 



2vr 



^/3 2cosh(2vrx/3) + 1 ' 



(19) 



/oo 
dxe lkx Ki(k) 
-00 
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with 

p -\k\ i , p -3fc/2-|fc|/2 1 _j_ „3fc/2-|fc|/2 

K (A;) = r, ^(fc) = ^—^ -r^, K 2 (k) = -rr, 

e k + 1 + e k e k + 1 + e k e k + 1 + e k 

\k\/2 — 3fe/2— |fc| 3fe/2— |fe| 

= fc— n ^(fe) = jrn i Zfc ,^s(fc) = fc— — _ fc , (20) 

e K + l + e K eHHe e e K + 1 + e K 

e -|fe|/2 + 2e lfc|/2 + e -3|fc|/2 
K6 ^) = fe , 1 , -fc ' 

e K + 1 + e K 
and the convolution is defined by 

J-OO ^ 

We introduce the rescaled eigenvalue A.'(v) = A(v)/[(ft+(v + i)<j)-{y — i)] which is useful, 

because of the simple asymptotic behaviour lim A (v) = const. At v = 0, the relation of 

v— »oo 

eigenvalue and rescaled eigenvalue is simply 

lnA(O) =mA'(0)-/3, 

i.e. it only amounts to a shift in the ground state energy. In terms of the auxiliary 
functions the rescaled eigenvalue reads 

lnA'(x) = -0e(x) + Vl * \nSi(x) + V x * hxS 2 {x) + V\ * \nS 3 (x) 
+ V 2 * In S 4 (x) + V 2 *ln S 5 (x) + V 2 *ln S 6 (x) 

with 

l + e~l fc l _ / vr 



From considering the low temperature limit by following Q, we conclude that there are 
two kinds of elementary excitations. One of them corresponds to the vector representation 
(si, §3,55), and the other to its conjugate (s 2 ,S4,se). Their bare energies are represented 
by ei,---e 6 . 
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Limiting cases 

Next, we take a suitable limit of the chemical potentials jXi, fi% and fi3, for which the 
spin-1/2 Heisenberg model is obtained. It is quite a natural result on the level of the 
Hamiltonian formalism, however less trivial on the level of the final NLIE, which supports 
the validity of our construction and the correctness of our calculations. Let us concentrate 
on the case of fj,± = // + hi, fj>3 = fj,' + /13 and \i2 = hi, where \h\\, \h 2 \, \h^\ <C |//|. First 
consider the limit fjf = 00. In this case, observing the constants in the driving terms of 
each auxiliary function, we conclude 

Sl = 0(l), S2 = 0(e-^'), s 3 = 0(1), (21) 
S4 = 0(e-^'), s 5 = 0(e^'), 8 6 = 0(e-M). (22) 



Therefore, we can regard 

S2 = S4 = s 6 = 0, 

s 5 ~ S 5 (23) 

that means the conjugate modes are suppressed. With the above ansatz, we can solve the 
equation for the ^-function as 

lnS 5 (x) = K 3 * (K - l)" 1 * {lnSi(x) + lnS 3 (x))+P(K - l)" 1 * V 2 {x) 



Substituting the last equation into (12) for s\ and S3 and performing the inverse Fourier 
transform, we get 

\np{x + i/2) = -2<!r(3$(x + i/2) + (3h/2 

/°° dv f°° dv 

JLk{x-y)\uP{y + i/2)- / ^-k(x - y + i) lnP(y - i/2), 
-00 ex J-00 ^vr 

\np{x-i/2) = 2irf3<5>(x-i/2) - (3h/2 (24) 

- f°° p.k(x-y-i)hxP(y + i/2)+ H p-k(x - y)lnP(y - i/2), 

J-00 ^vr J-00 ^7T 

. , n . % f°° , \nP{x + i/2) i f°° , lnP(x-i/2) 
lnA(x) = 2/? In 2 + (3[i+ - / dx K / dx 



2 J-00 sinh7r(x + i/2) 2 „/_oo sinh7r(x — i/2) 
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with 

p(x + i/2) = sx(x), p(x - i/2) = s 3 (x), P(v) = 1 + p(v), P(v) = 1 + p(v), 
^ , i 1 v fdx e ikx h 1 + h 2 + h 3 

* {v) = "^hW k{x) = J 2^TT^' h = hl - = 3 ' (25) 

where we have dropped the contribution of fi' to the potential [i. This NLIE is nothing 
but that for the spin 1/2-Heisenberg model ||. Thus, we see that the spin 1/2-Heisenberg 
model is obtained as a limit of the (3, 0)-US model directly in the two sets of NLIE's. 

Now, let us consider the opposite limit y! = — oo. In this case, we find the following 
simplifications 

s\ = s 3 = s 5 = 0, 

s 2 ~ S 2 . (26) 

Therefore, we obtain 

lnS 2 (x) = K 3 * {K - I)" 1 * (ln5 4 (x) + In 5 6 (s)) + P(K - l)" 1 * V x {x). (27) 



However, substituting (27) into (|T^), we see that 54, sq and A are constants. This fact is 
reasonable, because in the present case only the 2nd state can survive with finite energy. 
Hence, all physical degrees of freedom are frozen out at finite temperature. 

4 Numerical Analysis of the NLIE 

As an application of the obtained NLIE (jl^), we show numerical results for some physical 
quantities. Let us consider the entropy and the specific heat: 

S — m, C = T( B 4). (28) 



^dTJ \dT 

To avoid numerical differentiations, we simultaneously have solved the NLIE for the deriva- 
tives using relations like Q 

d wi , \ Si 8 1 
— In 1 + Si) = — -w^lnsi. 

Op 1 + Si op 

We have calculated the entropy and specific heat of the (3, 0)-US model with /xi = \x 2 = 
/i3 = numerically as shown in Fig. 2. 

II 



We can compare the result with that from the SL(3)i conformal field theory (CFT). 
From CFT and its finite size analysis, the free energy is predicted to be 

f = fo~^ + .... (29) 

Putting c = 2 and v = 2ir/3 [|0|, the low-temperature asymptotics is 

C = T. 

With a glance at Fig. 2, we see that the low temperature behaviour of the numerically 
determined specific heat has slope 1 as predicted by CFT. In a similar manner, we can 
calculate the specific heat in the presence of the chemical potential. First, we show results 
for the case fi\ = us = /u, fi2 = in Figs. 3a and 3b. As we have analyzed before, the 
specific heat in this case approaches that of the SU(2) spin-1/2 Heisenberg chain with zero 
magnetic field as [i —* 00. For the Heisenberg chain, we should put c = 1 and v = ir and 
get the asymptotics C = T/3. Indeed, the curves in Fig.3a show the expected tendency. 
Results for the specific heat for [i\ = ^3 = \i < 0, ^2 = are shown in Fig. 3b. Here we 
see a suppression of the specific heat data at fixed temperature for n — > — 00. This is the 
manifestation of the freezing of the system. Second, the case of [i\ = — [13 = h, 1x2 = is 
considered in Fig. 4. In this case, the rich sl(3) structure is most clearly exposed. 

At h = 4 in Fig. 4, we observe two structures (one shoulder and one maximum) of the 
specific heat owing to the fundamental mode and its conjugate. 

Let us define the magnetization M and the magnetic susceptibility x by 

d 2 f 



M = --f, x- 



h=0 



dh 2 



h=0 



dh' dh 

respectively. The results are shown in Figs. 5a, b and 6. From Figs. 5 we read off the 
existence of two critical field h c % = 0.94 and h C 2 = 4. For h > h C 2 only the a = 1 state can 
survive at T = 0, thus the groundstate magnetization is maximal, i.e. M = 1. This is in 
agreement with analytic considerations. The lower field h c \ is the value, above which only 
the q = 1 and a = 2 states can exist. Note that for these critical fields the magnetization 
M(T) shows a square root behaviour at low T. 
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The magnetic susceptibility is given in Fig. 6. The curve is qualitatively similar to that 
of SU(2) Heisenberg model however with x(0) = 3/V 2 . The susceptibilities at the lowest 
temperatures shown in Fig. 6 are still about 10% above the groundstate value. The origin 
of this singular behaviour of the susceptibility at T = are 1/logT corrections similar to 
1 17, and will be discussed elsewhere. 



5 Summary 

Before closing this paper, we would like to mention another application and generaliza- 
tions of the NLIE. First, we can use our NLIE in order to calculate the characters of the 
SL(3)i Kac- Moody algebra. Let us remember that our NLIE contains a natural notion of 



spinons and chemical potentials |19|| -[21|. This fact enables us to present the characters 



in terms of the spinons and chemical potentials [22]. Details will be published elsewhere. 
Second, up to our knowledge the complete strings in the anisotropic sl(n) (n > 2) models 
with trigonometric and elliptic i?-matrices have not been constructed yet. However, one 
of the present authors has shown that the QTM method is applicable not only to the 
isotopic XXX Heisenberg chain but also to the anisotropic XX Z and XYZ versions [||. 
Therefore, we expect that a combination of the anti-symmetric fusion procedure and the 
NLIE will be generalizable to the anisotropic s/(3) models. 
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Figure Captions 

Fig.l i?-matrix for the PS model. 

Fig. 2 Entropy and specific heat in the absence of the chemical potentials. 

Fig. 3a Specific heat for fii = ^3 = fi > 0. 

Fig. 3b Specific heat for fii = ^3 = fi < 0. 

Fig. 4 Specific heat with fii = —fi 3 = h > 0. 

Fig. 5a Magnetization for fii = —^3 = h C 2, ^2 = 0. 

Fig. 5b Magnetization for fii = —^3 = h c i, fi2 = 0. 

Fig. 6 Magnetic susceptibility. 
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